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tion, is found. The vacuum setting is done by the 6D Higgs potential and 
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from 6D to 4D by use of the soliton solution (brane world) . It is regarded as 
a higher dimensional version of the Randall-Sundrum 5D model. The mass 
hierarchy is examined. Especially the geometrical see-saw mass relation is 
obtained. Some physical quantities in 4D world such as the Planck mass, 
the cosmological constant, and fermion masses are focussed. Comparison 
with the 5D model is made. 
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1 Introduction 



The higher dimensional approach is a natural way to analyze the 4D physics 
in the geometrical standpoint. The history traces back to Kaluza-Klein in 
1921 Q, g where the electro -magnetic force is amalgamated into the gravita- 
tional force in the 5D theory. Stimulated by the recent development of the 
string and D-brane theories, a new type compactification mechanism was 
invented by Randall and Sundrum^, ||]. The domain wall configuration in 
5D space-time, which is a kink solution in the extra dimension, is exploited. 
Before Randall-Sundrum(RS), the supersymmetry is the popular way to ex- 
tend the standard model. The D-brane inspired model has provided us with 
new possibilities for the extension. It has some advantages in the hierarchy 
problem and the chiral problem. 

A moderate view to the string or D-brane theory is to regard the for- 
malism as a regularization of the field theory. In this standpoint, the ex- 
tendedness, in the space-time, of an object is the essence of the singular-free 
behavior. The non-commutative geometry, which is a recent trend in the 
field theory, can be regarded as a systematic way to introduce the extend- 
edness. In the ordinary (commutative) geometry, the soliton configuration 
can do the same role. We present a 6D soliton solution, and show that it 
provides a new dimensional reduction mechanism. 

In the research of the domain world physics, the main theoretical fo- 
cuses, at present, are 1) Consistency, 2) Stability and 3) Localization. The 
first one is the basic thing in all theoretical analysis. In the analysis of the 
domain world physics, the delta-function type special configuration (singu- 
lar space) is often assumed without properly solving the field equation with 
some boundary condition. This is a main origin of the necessity of the con- 
sistency check. In some limitted situation, it can be accepted and is indeed 
useful in that the whole physical image can be easily introduced. The ap- 
proach is, however, not admitted as the right configuration in the general 
case. It is sometimes discussed that some no-go theorem exists in some 
models [||, ||, 0. Especially the rigorous treatment is required in the analysis 
related to the cosmological constant and the boundary condition ||, ^, [TO]. 
Recently in order to check such consistency for a wide-range brane-world 
models, some useful sum rule is presented[ll|. The present approach treats 



the configuration in the general setting and the model is based on the field 
theory. Therefore it basically passes the consistency check. Next item is the 
stability problem. The solution should satisfy the stability condition. It is 
generally rather difficult to analyze because the solution should be examined 
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from some general configuration depending on parameters such as the dis- 
tance between the branes. One interesting idea is proposed by Goldberger 
and Wise|12[], where the radion field is utilized. In the present approach, 
the stability is guaranteed by the boundary condition. This is the popular 
situation in the soliton physics. The third focus is whether the ordinary 
fields (matter fields, gauge fields, graviton, etc.) are localized in some nar- 
row region in the extra space. This is necessary to be consistent with the 
fact that we have so far not observed the extra dimensions (requirement for 
the dimensional reduction). In the present analysis, no real matter fields are 
taken into account. Only 6D Higgs scalar is introduced to define the (classi- 
cal) vacuum clearly. The localization of other fields, in the present analysis, 
is assumed valid based on other references [13, 14] and on the analogy to the 
lattice fermion situation. 

We have some reasons to investigate the higher dimensional general- 
ization of the RS-model. Firstly some references suggest the difficulty of 
making a realistic stable 5D domain world ]|, |(|. One of main causes looks 
its odd dimensionality. Secondly, in the solitonic approach, various extended 
objects (wall, string, monopole, ...) appear[15]. It is a natural direction to 
investigate the possibility of other type configuration than wall. Thirdly in 
the work by Callan and Harvey Jl^], they interpret the anomaly phenomena 
in 4D world in the higher dimensions. They exploit the physics of fermion 
zero modes on strings in 2n+2 dimensions and domain walls in 2n+l dimen- 
sions. The result is consistent with the anomaly relations between different 
dimensions by Alvarez Gaume and Witten[17|. 6D model is necessary, in 
addition to 5D, to understand the whole anomaly structure. We consider 
the brane world model has great advantages about the anomaly and chi- 
ral problems. Finally, in the present rapid progress of the noncommutative 
geometry, the extra space has even dimensions. We suppose the present 
approach has something common with such approach. 

The 6D models of the hierarchy problem have already been discussed in 
various waysfl8|, |l9| , p0| , f2~l| . Main different points from those works are 1) 
the present approach puts emphasis on the thickness of the configuration, 
2) the 6D(bulk) Higgs fields are introduced to clearly define the (classical) 
vacua, which are necessary for the soliton boundary conditions, 3) we treat 
the warp factors and the scalar (Higgs) fields on the equal footing, 4) the 
boundary condition is systematically solved and 5) supersymmetry is not 
taken into account. 

In Sec. 2 we introduce the 6D model and derive a solution of the 6D 
classical Einstein equation in Sec. 3. It is a one parameter family of soliton 
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solution with the pole configuration. The asymptotic forms, in the dimen- 
sional reduction, of the vacuum parameters are obtained in Sec. 4. In Sec. 5, 
some physical quantities, such as the 4D Planck mass , the 4D cosmological 
constant and fermion masses, are explained. The see-saw relation is obtained 
there. Some order estimation is also done. The precise form of the solution 
is obtained by analytically solving the boundary condition in Sec. 6. It is 
demonstrated that the vacuum parameters are fixed to be some numbers by 
one input (boundary) condition. The obtained solution is further confirmed 
by the numerical method in Sec. 7. Finally we conclude in Sec. 8. In app.A 
we explain, in detail, the derivation of the solution treated in Sec. 3, . Some 
detail about the numerical analysis of Sec. 7 is explained in App.B. 

2 Model set-up 

We start with the 6D gravitational theory, where the metric is Lorentzian, 
with the 6D Higgs potential. 

S[G A b, $] = J d 6 xV^G(-^M 4 R - G AB d A <f>*d B $ - V($*, $)) , 

V(^^) = ^\ 2 - Vq 2 ) 2 + A , (1) 

where X (A = 0,1,2,3,4,5) is the 6D coordinates and we also use the 
notation (X A ) = (x^, p, <p), p = 0,1,2,3. x M 's are regarded as our world 
coordinates, whereas (X 4 ,X 5 ) = (p,<p) the extra ones. The extra (spacial) 
coordinates are taken to be polar ones: < p < oo, < ip < 2tt. p will be 
regarded as a "radius" and if as an "angle" around the 4D pole- world. See 
Fig.l. The signature of the 6D metric Gab is ( — h + + ++)• $ is the 6D 
complex scalar field, G = det Gab, R is the 6D Riemannian scalar curvature. 
M(> 0) is the 6D Planck mass and is regarded as the fundamental scale of 
this dimensional reduction scenario. V ($*,$) is the Higgs potential and 
defines the (classical) vacuum in the 6D world. The three parameters A, vq 
and A in the potential V are called here vacuum parameters. A(> 0) is 
a coupling, vq(> 0) is the Higgs field vacuum expectation value, and A is 
the 6D cosmological constant. It is later shown that the sign of A must 
be negative for the proposed pole vacuum configuration. We take the line 
element shown below. 3 

ds 2 = Q-^n^dx^dx" + dp 2 + p 2 e~ 2 ^ p) dip 2 , 
3 The same type metric was also taken in Ref.M, p0[. 
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< p < oo , < if < 2vr , (2) 

where = diag(— 1, 1, 1, 1). This is a natural 6D minimal-extension of 
the original 5D model by Randall- Sundrum ||. For a = lu = 0, ds 2 is 
the 6D Minkowski (flat) space. In this choice, the 4D Poincare invariance 
is preserved. Two "warp" factors e~ 2a ^ and e~ 2u) ^ appear and play an 
important role throughout this paper. Note that, for the fixed p case {dp = 
0), the metric is the Weyl transformation of the product space of the flat 
(Minkowski) space, rj^dx^dx" , and the circle S , p 2 dip 2 . The two Weyl 
factors serve for scaling each space. The coordinate p is regarded as the 
scaling parameter. 

From (H) we can read 

\ / e" 2 ^ \ 
Gab \= 1 . (3) 

] \ p 2 e- 2 " J 

The 6D Riemannian scalar curvature R is given by 

R = 4(-2a" + 5{a') 2 - 2—) + 8wV + 2(-u" + {to') 2 - 2—) , (4) 

P P 

where ' = ■$-• We will later require there are no curvature singularities any- 
where because the configuration considered has the finite thickness (radius 
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of the pole) which should smooth the curvature. Since M-dependence can 
be absorbed by a simple scaling ($ = M 2< S,v = M 2 vo,X = M~ 2 X,A = 
M 6 A, = M~ l x^, p = M~ l p, (p = <p), it is sometimes useful to take 

M = 1 . (5) 

We will, however, explicitly write M as much as possible in order for the 
reader not to bother about the superficial dimensional mismatch. 



3 A solution 

Let us solve the 6D Einstein equation. 



M\R M n - \GmnR) = 
-d M <S>* d N $ - d N $* d M <S> + G M N{G KL d K ^* d L <& + V) , 

V 2 * = l^ • (6) 

As for the complex scalar field $, first of all, it should satisfy the periodicity 
condition: $(x, p,(p) = <I>(x, p, <p + 2tt). Following Callan and Harvey [pi, we 
consider a simple case that <3? depends only on the extra coordinates ip and 
p as 

$ m (p,<p) =P{py mLp , m = 0,±l,±2,--- . (7) 
Then the Einstein equation @ reduces to 
M = p, N = v 

3cJ " _ 6(T ' 2 + 3- - 3a'u' + J' - [Jf + 2- 
P P 



2lo 

M-\P' 2 +m 2 ^P 2 + V) , (8) 
P l 



M = p,N = p 



_ 6cJ ' 2 + 4^ = M - A {-P' 2 +m 2 -^P 2 + V) , (9) 

P P 

M = <p,N = <p 

Eq.(C) 4a" - Wa' 2 = M~\P' 2 - m 2 -^P 2 + V) , (10) 



Matter Eq. 



p 2u) 1 \ 

P " _ m 2 —P - Ua' - - + J)P' = -(P 2 - v 2 )P . (11) 
P 2 P 2 
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The matter equation (O) is derived from other three equations. From 
(||) and @, and from (||) and ([k]), we obtain two convenient equations as 
follows. 



Eq.(A) Za" -- + a'J + u" -(Jf + 2- = 2M- i P' 2 , (12) 
P P 
i 

Eq.(B) -16a' 2 +4 Aa'J + 4a" = 2M~W . (13) 

P 

The three equations Eq.(A),Eq.(B) and Eq.(C) may be considered as those 
to be solved. They have three unknown functions P(p),a(p),u(p) and one 
unknown integer parameter m. We call this set of the three equations reduced 
Einstein equations (REE). 

Now we must consider the boundary condition(asymptotic behavior), at 
p = oo(infrared region), for P(p),a(p) and oj{p). As for P(p), we naturally 
take 

p^oo , P(p)^±v . (14) 

The plural sign ± comes from the fact that the Higgs potential V(&) is the 
even function of <£. Later soon, we will choose the plus sign above by use of 
the freedom of an arbitrary parameter eo( defined later ). As for other two, 
we assume (from the "experience" in 5D Randall-Sundrum model||, ^j, [|) 
a' — > a(const), u/ — > 6(const) as p — > oo. Then from the REE, we can deduce 




As p — > oo , a' — > a , u>' — > a , a = ±y^-M ^ . (15) 

As for the plural sign above, we take the plus one which damps the scaling 
behavior (^) in the asymptotic region. Note that the present asymptotic 
requirement demands the isotropic property around the p = axis (m = 0), 
that is, $ should not depend on the direction ip, that is, the pole configu- 
ration. See Fig.l. We take the plus sign above in order for this mechanism 
to work, in the infrared region, for exponentially suppressing the 4D part of 
the metric (g). In the above result we must have the condition: 

A < . (16) 

4 It reflects the following well-known fact. The matter equation V 2< 3? = j^r in (^) can 
also be derived from the (6D) energy-momentum conservation M 4 V iu (Rmn — \GmnR) = 
V M T M N = 0, where T M n =RHS of the first equation of (ph. 
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The geometry should be Anti de Sitter in the (infrared) asymptotic region. 
We should further note that, in the infrared region above, the 6D scalar 
curvature is a positive constant: R — > 30a 2 = —3AM -4 > 0. 

Next we must consider the boundary condition at p = +0(ultra-violet 
region). Assuming the following things: 1) As p — > +0, the three unknown 
functions a',u' and P goes like a' — > sp a , u/ — > wp b , P — > xp c where 
s,w and x are non-zero constants (s / 0,w / 0,i / 0); 2) The 6D scalar 
curvature R is regular everywhere; 3)^| a = b, we can deduce, from REE, 
that there are two cases: 

Case 1 a = b = 1 , c = , s = — (x 2 - v 2 ) 2 + - (17) 

16 4 

Case 2 o = 6 = l,c>l,s = — v 4 + t ( 18 ) 

16 4 

If we put off the condition x ^ in fll7|) and take x = 0, eq.(|l7|) reduces to 
&■ 

Let us take the following form for a' (p) , ui' (p) and P(p) as a solution ^| . 

oo 

a'(p) = a y , C2n+ \, {tanh(M} 2w+1 , 
W ^ (2n + l)! 1 1 In 

OO T 

= a Y , d2 " + \, {tanh(M} 2n+1 , 
W ^ (2n + l)! 1 1 

oo 

= E 7Sj{ tanh (M} 2n , (19) 

n=0 ^ ' 

where a = +J ^-M~ 2 and c's, d's and e's are coefficient-constants (inde- 
pendent of p) to be determined. Because REE have the "formal" symmetry 
of parity: p — > —p, we know the solutions are parity odd or even functions. 
(Note that < p < oo, @. ) In fli~9|), a' and u>' are composed of odd 
powers of tanh(fcp), whereas P is of even powers. The choice is deduced by 
the behavior at the ultra-violet region dHH). The expansion (|T9| ) could 



be regarded as a sort of " Kaluza-Klein mode expansion" . f\ In the present 

5 The assumption a = 6 is based on the expected " duality" between IR and UV 
behaviors. Note that, in the IR, a' and lj' behave similarly in (p_5|). 

6 Some properies of the form of the series (|l^) in relation to the important equation 
( |4l| ) is described in App.A. 

7 In comparison with the work by Goldberger and Wise|%|, not only the scalar field 
P but also the warp factors a' and u)' are "Kaluza-Klein expanded". 
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case, however, there is no periodicity in the p-coordinate. We neither have 
the translation invariance p — > p + const., which should be compared with 
5D case[Q, ||. A new mass scale k(> 0) is introduced here to make the 
quantity kp dimensionless. The physical meaning of 1/k is the "thickness" 
of the pole. The parameter k, with M and p c (defined later), plays a central 
role in this dimensional reduction scenario. We call them, (k,M,p c ), funda- 
mental parameters. The distortion of 6D space-time by the existence of the 
pole should be sufficiently small so that the quantum effect of 6D gravity 
can be ignored and the present classical analysis is valid. This requires the 
condition |EJ 

k < M . (20) 



The boundary conditions (|TjJT^) with the + choice in the plural signs, 
require the coefficient-constants to have the following constraints 



OO OO 7 oo 

.. _ \ ^ C2n+1 1 _ \ - «2n+l -. _ \ - &2n 



(21) 



, which are obtained by considering the asymptotic behaviors p — > oo in (]T9) . 
We will use these constraints in Sec. 6. In App.A, we explain the determina- 
tion of all coefficients in ( [T9D from the REE. All ones are fixed except one 
free parameter eo(= P(0)/vq). ^ The first two orders are concretely given 



as 



<* = w{lW(l-eo 2 ) 2 + A} 
di = -fci , 
eo : free parameter, 

c 3 = I ^M~W(1 ~ eo 2 ) 2 + ci(2 + 5f d) 
d 3 = -|ci(l + 5fci) , 



3^1^ i 

e 2 = -|^e (l-eo 2 ) , 

(22) 



Here we note that those parts of the UV boundary condition (|17|jl8| ), which 
are not taken into account in the starting assumption (|l^), are satisfied by 

8 The corresponding global symmetry looks some 'deformation' of the constant transfor- 
mation: P(p) — > P(p) + eo- The symmetry is valid for the special case: m = 0, A = 0. This 
situation should be compared with the 5D case where the constant translation symmetry 
for a coordinate exists^ ^j. 
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the above solution: [Case 1] eo ^ 0,x = t>oeo,s = ac\\ [Case2] eo = 0, s = 
ac\ . The third order is given by 



e 4 xM- { ^e 2 + l^e (l-eo 2 )} 
a 12 ka 

4 a 5 

= (-7TC3 + 2ci) + - • -(dic 3 + d 3 ci + 2cic 3 ) 
9 k 3 

Av 4 l, /f -4/ 1 o 2n 2 lR W 2 ,,,-4 2 

— ; — -M (1 — 3e )e 2 +6 M e 2 , 

fca 4 a 

5 Aw 4 w 4 ,„ 2N 100 

c 5 = ~J^-^ M e o(! " eo 2 )e 4 + — c 3 + 4ci 

a 10,, , . 5 Afn 4 , , a, 9, 9 

+ i ' y ( 103 + 3Cl + lCs) " ii^r ( " )e2 ' 

120 r 2 Aw 2 13 1 11 , 2 , 

d 5 = -2c 5 + — {-^ L M- 4 e 2 e 4 + — c 3 - -ci + — d 3 + -d x 
7 3a 9 5 18 5 

"I ' ^ lC3 + d3Cl " 2dlds) ~ A ^~ M ~ 4e22} (23) 

Note that, in the first equation, e 4 is written by lower-order terms 
(eo, e 2 ; ci, c 3 ; d\, d 3 ). Putting this result of e 4 into the second equation, we 
see C5 is written by the lower-order terms. Similarly to the third equation 
with respect to d§. 

The general terms (c 2n+ i, d 2n +i, e 2n ), n > 3 are obtained in App.A. They 
are expressed by the coupled linear equations fl57| ) with respect to them and 
can be solved recursively in the increasing order of n. All coefficients are 
expressed by four parameters A, vo,A and eo- (We may take M = 1 as 
explained before. We may also take k = 1, keeping the relation (20) in 



mind. This is because the thickness parameter k can be absorbed by the 
change : c 2n +i -> fc 2 „+i,d 2n+ i -> |d 2n +i,A -> ^M _4 ,A -> ^M 4 . We 
will, however, keep fc as much as possible in order for the reader not to bother 
about the superficial dimensional mismatch.) In the above results, we see 
&2n{n > 1) is the odd function of eo, whereas c 2n+ i and d 2n+ i are the even 
functions. Under the sign change eo — > — eo, P changes the sign, whereas a' 
and {J do not. As commented before, we can choose the plus sign in ( [14]) 
by use of this freedom. The above four parameters have three constraints 
( pl| ) from the boundary condition at the infrared infinity. Hence the present 
solution is one-parameter family solution. We will solve the constraints in 
Sec. 6. As a final note in this section, we point out that, in comparison with 
5D RS-model[||, 0], there appears no lower-bound for A. 
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4 Meaning of the Solution 



Let us examine the obtained solution in relation to the two general theo- 
rems: Goldstone's theorem and the Derrick's theorem. 



(i) Spontaneous Symmetry Breakdown 

The present model (|]) has the discrete symmetry: 

$ -> , -> -$* , (24) 

and the continuous (global) symmetry: 

$ _^ e i0 $ , -► e _i0 $* , (25) 



where 8 is a constant parameter. (|24| ) is a special case of (25) and appears 



also in the case of the real scalar field. If the vacuum breaks the continuous 



symmetry, a massless scalar boson appears (Goldstone's theorem [23, |24|j). 
In the present case, however, we confine the complex scalar field <E> to the real 
field by taking only the m = mode in (^). (In fact we can start from the 
real scalar in (Q) and obtain the same result.) This is required from the 
demand that the IR boundary condition of the warp factors is the same as 
the 5D RS-model. Hence we should understand that, in the present solution, 
the discrete symmetry (p3|) is spontaneously broken, not the continuous one 
(p5|). No Goldstone boson appears, and the situation is similar to the the 
domain wall case of 5D model g |. When we want to consider the breaking 
of the continuous symmetry, the standard way is to take the Abelian Higgs 
model[^5|] where the Goldstone boson is 'eaten' f\ by the gauge field (Higgs 



phenomenon). Other approach is also examined in [26|. 



(ii)Role of the Warp Factors 

After taking the 'warped' metric (|2|) and assuming the form of the complex 
scalar field (^) with m = 0, the action reduces to 

S = 2vr / d 4 x dp pe- 4c7 - w {e 2a <S>*<i> - P' 2 + 2(2a" - 5a' 2 + 2—) 
J P 

-Au'a' + (u" - J 2 + 2—) - -(P 2 - v 2 ) 2 - A} (26) 

P 4 

where 6 = 4r which vanishes in the present static solution. If we omit 
the kinetic term of in the above integral formula, the remaing part is 
9 The massless Goldstone mode is replaced by one mode of the massive gauge field. 



11 



regarded as (-l)x (Energy of the static solution). Its finiteness is apparent 
from the factor of e _4<J_a; and the regular behavior of P, a' and to' near p = 0. 
The expression (|26|) looks like an interacting theory of scalars. As for the 
localized solutions in scalar theories, there exists the well-known theorem 
called Derrick's theorem^, p3| . It says, in the case of the space dimensions 
higher than or equal to 2, the only non-singular time-independent solutions 
of finite energy are the ground states. The present case escapes this trivial 
situation in the following way. As shown in fl26| ) , the interraction forms of the 
fields a(p),u>(p) are not those of scalars but those of dilatons. Especially the 
linear terms (a", a' / p, lo", uj 1 / p) appear. They break the positivity property 
used in the proof of the Derrick's theorem. 



5 Vacuum parameters: M and /c-dependence in the 
dimensional reduction 

Let us examine the behavior of the vacuum parameters (A, vq,,X) near the 
4D world limit (thin pole limit): k — > oo(the dimensional reduction). This 
should be taken consistently with (|2C). We will specify the above limit in 
the more well-defined way later. In (19)), we note that {tanh(/c/j)} 2n+1 —* 
0(kp) , {tanh(/cp)} 2n — > 6{kp) , p > as k — > oo . Assuming the convergence 
of the infinite series, we can conclude that c\ ~ 0(k°) , e2 ~ 0(k°) as 
k — > oo. (eo is a free parameter.) From the former relation, using the explicit 
result of ci (|22"|), we know —A ~ A^o 4 ~ M A ka. Because a ~ y/— AM -2 , 
this says V~ A ~ kM 2 and Xvo A ~ k 2 M A . From the latter one, using the 
explicit result of &2 (^2[), we know At>o 2 ~ k 2 . Hence we obtain 



- A ~ M 4 k 2 , vq ~ M 2 , A ~ M~ 4 k 2 as k -> oo . (27) 

These are the leading behavior of the vacuum parameters in the dimensional 
reduction from 6D to 4D as k — > oo. As for the fc-dependence, the above 
result is the same as the 5D model of Randall and Sundrum|3l Hi. The more 



precise forms of (27) will be obtained, in Sec. 6, using the constraints (21). 



6 Physical constants and See-Saw relation 

In order to express some physical scales in terms of the fundamental param- 
eters M, k and p c ( to be introduced soon), we consider the case that the 4D 
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geometry is slightly fluctuating around the Minkowski (flat) space. 
ds 2 = e-^g^dx^dx" + dp 2 + pV 2aj(p) V , 

9nv = V^u + h^v , h pv ~ O(-) . (28) 
The leading order O(k ) results of the previous section remain valid. 

6.1 The Planck mass 

The gravitational part of 6D action (|l]) reduces to 4D action as 

J d 6 xV^GM 4 R ~ M 4 dip J* dp pe -Mp)^(p) J d A x^=g~R + ■■■ , (29) 

where the infrared regularization parameter p c is introduced. p c specifies the 
size of the extra 2D space. Using the asymptotic forms, a ~ ap, to ~ ap as 
p — > oo and a = J ^M~ 2 ~ k as — ► oo, we can evaluate the order of M p i 
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as0 



MJ~MHv\ d ppe- 3a P~ ~ , (30) 



l p l ~ J Q a Pf Je ~ ~ju2~ 

where we have used the reduction condition: 

kp c > 1 . (31) 



2 



The result ( |30|) is different from 5D model of Randall-Sunsrum [|3j, |8j]: 

M 3 

The above condition fl3l|) should be interpreted as the precise (well reg- 
ularized) definition of k — > oo (dimensional reduction) used so far. We 
note p c dependence in ( |30"1) is negligible for /cp c 3> 1. This behavior shows 
the distinguished contrast with the Kaluza-Klein reduction of 6D to 4D 
(M p i 2 ~ M 4 (p c ) 2 ) as stressed in |j for the 5D case. 

Writing @ as 

-W~T ' (32) 

we notice this mass relation is the geometrical see-saw relation corresponding 
to the matrix : 



M 

M M p i 



(33) 



A similar relation is derived in other 6D model 1201 
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This provides the geometrical approach to the see-saw mechanism |2^, ^] 
which is usually explained by the diagonalization of the (neutrino) mass 
matrix. ( See a textbook|3^].) For some discussion about the neutrino 
mass, using the above relation, see the final section. 

6.2 The cosmological term 

The cosmological part of (||) reduces to 4D action as 

J (fxV^GA ~ A dip dp P e- 4a ^-"W J d 4 x^ = A 4d J d A x^ 



A 

7? 



A id ~ ^ ~ —M A < , kr c > 1 



A^ is the cosmological term in the 4D space-time. It does not, like M p [, 
depend on p c strongly. The above result differs from RS-model in that k- 
dependence disappears, (cf. A 4 ^ ~ —M 3 k in RS-model ||.) The result says 
the 4D space-time should also be anti de Sitter. 

6.3 Numerical fitting of M, k and p c 

Let us examine what orders of values should we take for the fundamental 
parameters M and k. ( p c is later fixed by the information of the 4D fermion 

parameter = A^/Mpi 2 [^|] has the relation: 



masses. ) Using the value M p i ~ 10 19 GeV , the "rescaled" cosmological 



A M ~ k ~ M 2 x 1(T 19 GeV , (35) 



where the relations (^) and ( |34| ) are used. The unit of mass is GeV and 
this unit is taken in the following. The observed value of A^ is not definite, 



even for its sign. [32] If we take into account the quantum effect, the value of 



A^ could run along the renormalization[33|. Furthermore if we consider the 
parameter A^ represents some "effective" value averaging over all matter 
fields, the value, no doubt, changes during the evolution of the universe. 
(Note the model (||) has no (ordinary) matter fields.) Therefore, instead of 
specifying A^, it is useful to consider various possible cases of A^ ~ —k 2 . 

Some typical cases are 1) (k = 1(T 41 ,M = 1CT 11 ), 2) (k = 1(T 13 ,M = 
10 3 ) 3) (k = 10, M = 10 10 ) 4) (k = 10 4 ,M = 10 1L5 ) and 5) (k = 10 19 ,M = 
10 19 ). Case 1) gives the most plausible present value of the cosmological 
constant. The wall thickness 1/k = 10 41 [GeV _1 ] , however, is the radius of 
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the present universe. This implies the extra dimensional effect appears at 
the cosmological scale, which should be abandoned in the usual sense. For 
further discussion, see the final section. Case 2) gives 1/k = 10 13 GeV -1 ~ 



1mm which is the minimum length at which the Newton's law is checked[34, 



35 1 . Usually k should be larger than this value so that we keep the observed 
Newton's law. 5) is an extreme case M = M p i. The fundamental scale M 
is given by the Planck mass. In this case, p c S> 1/k = 1/M p i is acceptable, 



whereas y — ~ M p i is completely inconsistent with the experiment and 
requires explanation. Most crucially the condition (20) breaks down. Cases 



3) and 4) are some intermediate cases which are acceptable except for the 
cosmological constant. They will be used in the next paragraph. At present 
any choice of (k, M) looks to have some trouble if we take into account 
the cosmological constant. We consider the observed cosmological constant 
(10 GeV) should be explained by some unknown mechanism. 

As in the Callan and Harvey's paper [ 16 1, we can have the 4D massless 
chiral fermion bound to the wall by introducing 6D Dirac fermion ip into 

(©• 

S[G A b,<S>}+ J d 6 xV Z G(^i> + g$(Re<5>-ir 7 lm<I>)4>) . (36) 

If we regulate the extra axis by the finite range < y < p c , the 4D fermion is 
expected to have a small mass mj ~ ke~ kpc (This phenomenon is known in 
the condensed matter physics as the surface mode. It is also confirmed in the 
two- walls configuration of the 5D lattice theory |J6|, |37|). If we take the case 

3) in Subsec.5.3 (k = 10, M = 10 10 ) and regard the 4D fermion as a neutrino 
{m v ~ 10 -11 - 10~ 9 GeV), we obtain p c = 2.76 - 2.30GeV _1 . If we take case 

4) (k = 10 4 ,M = 10 115 ), we obtain p c = (3.45 - 2.99) x lO^GeV -1 . 
When the quarks or other leptons (m q ,mi ~ 10~ 3 — 10 2 GeV) are taken 
as the 4D fermion, and take the case 4) in Subsec.5.3, we obtain p c — 
(1.61 — 0.461) x 10 _3 GeV _1 . It is quite a fascinating idea to identify the 
chiral fermion zero mode bound to the pole with the neutrinos, quarks or 
other leptons. 

6.4 The valid region of the present approach 

The condition on k in the present model, from ( f20|) and (|3l|), is given as 



— < k < M . (37) 

Pc 
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In ref.||, it is pointed out that the present mechanism has the similarity to 
the chiral fermion determinant. The above relation corresponds to \k^\ <C 
Mp <C i in ref.[p8|, |39| , Especially the thickness parameter k corresponds 
to 5D fermion mass parameter Mp which is a very important parameter, 
in the lattice simulation, to give good numerical results for the soft-hadron 
physics (such as the pion mass). 



7 Precise form of Solution (Analytical Approach) 

Let us determine the precise form of the solution by solving the boundary 
condition (^). An interesting aspect of the present solution is that some 
family of vacua is selected as the consistent (classical) configuration. In 
other words, the vacuum parameters are fixed to be some numbers (if one 
free parameter is fixed). Let us determine the precise form of (|27| ) and (|i~9l ) 
using the three constraints (21). The coefficients (c2 n +i, ^n+l, e-2n) of ( |l9| ) 
are solved, in (p2]) and (23), in terms of four parameters A,i>o,A and eo. 
Therefore the present solution is one parameter family solution. The precise 
forms are obtained by the p-expansion for the case kp c > 1 as 



-A = M 2 k(a + + •••) = M 2 ^ 



( k Pc) 2 ^ (kpc 

= M*k\8o + + •••) = M 4 k 2 f) *" 



\2n 



( k Pc) 2 {kpc 



,2n 



where a's,S's and (3's are some numerical (real) numbers to be consistently 
chosen using (|2l|). 

If we assume the infinite series of (|2l]) converge sufficiently rapidly,we 
can safely truncate it at the first few terms. In order to demonstrate how 
the vacuum parameters are fixed, we take into account up to n = 2 in (]2l| ) 
and the leading order in (|38|). 11 We present two sample solutions: 



Solution 1 eo = — 0.5(input), 
(a ,A),<5o) = (2.507019,1.113947,29.14319) ; 



11 In ref.[y], the improved calculation results are obtained. It takes into account terms 
up to n=6th order. 
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Solution 2 eo = — 0.8(input), 
(a ,A),<5o) = (1.864786,0.6424929,15.90691). 



(39) 

For each solution, the expansion coefficients are obtained as 
Solution 1 

(ci,cs,cb) = (-0.689617,13.1849,-60.9444), 
(di,d 3 ,d 5 ) = (0.459745,-1.59403,96.7112), 
(e ,e2,e 4 ) = (-0.5,2.45270,6.56764); 

Solution 2 

( Cl)C3)C5 ) = (-1.25575,7.33112,124.067), 
(d u d 3 ,d 5 ) = (0.837165,-4.52497,110.039), 
(e , e2 ,e 4 ) = (-0.8,1.78258,21.809). 

(40) 

For the solution 2, we plot P(p),a'(p),uj'(p) and R{p) in Fig.2, Fig.3, Fig.4 
and Fig. 5, respectively. We stress that the Riemann scalar curvature R is ev- 
erywhere non-singular. Note that R/a 2 — > —12^/lOcx/A, A = \/—A/(M 2 k) 
as p — > +0 (ultraviolet region). (Purely for the technical reason, Fig. 5 do 
not show the value of R at p = 0.) As for the infrared region, we know 
-R/a 2 — ► 30 as p — ► oo. 

Note the following points: 1) The above results are valid for general 
fundamental parameters (k,M,p c ) except the thin-pole condition kp c ~S> 1; 
2) For each given value of eo, the solution looks unique (when a solution 
exists) as far as the reasonable range of (ao, A)>^o) is concerned; 3) For the 
choice of the plus boundary for the Higgs field P (the plus sign in Ql4])), we 
do not find solutions for positive eo in the reasonable range of parameters; 4) 
For one given input value ( eo in the above case. The 6D Higgs mass (in unit 
k) \/\vo' 2 /k ~ y/do/flo is another example of such value.) the cosmological 
constant (in unit M 2 k) \/—A/M 2 k ~ ao is fixed ; 5) The above very 

12 When we see the convergence behavior of the series appearing in (jis]) and (^), we 
should take into account the factorial factors which divide the all coefficients c's, d's and 
e's. 

13 In some references, this is mistakenly called "fine-tuning". The cosmological constant 
is here not fine-tuned but fixed from the (classical) field equation and the input value. 
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Fig.2 The Higgs Field P(p)/v , ©■ Solution 2 in @ and @. 

Horizontal axis: kp. 



fine results are obtained by the standard method of the numerical calculus: 
Newton method. 

Our solution has one free parameter. Using this freedom we can ad- 
just one of the three vacuum parameters in the way the observed physical 
values are explained. In (|39|), we take eo = P(0)/vq as the adjustment. 
Solution 2 has the vacuum expectation value vqM~ 2 ~ yffio = 0.802, the 
cosmological constant AM" 4 !;" 2 ~ — «o 2 = —3.48 and the 6D Higgs mass 
V \vQ 2 k~ l ~ \/Sq/(3o = 4.98. We notice the dimensionless vacuum parame- 
ters, vqM~ 2 , Ak~ 2 M~ 4 and Xk~ 2 M , are specified only by the value of kp c 
and one input data, say, eo- When we further specify k and M, as consid- 
ered in Subsec 5.3, the values Vq, A and A are themselves obtained. Any 
higher-order, in principle, can be obtained by the ^-expansion. 



8 Numerical Solution 

Instead of the previous analytic approach based on the assumed function 
form (|T9|), we can directly solve the coupled differential equations Eq.(A), 
Eq.(B) and Eq.(C) ((12), (13) and (|i~0|)) in the numerical method (Runge- 



Kutta method). The difficult point of this approach is the right choice of 
the vacuum parameters. The three parameters should be correctly chosen 
corresponding to a given boundary condition. (In the previous section, say, 
the condition is the initial value of P: P(p = 0)/vq = eo-) Another difficulty 
comes from the technical reason in the numerical analysis. In the coupled 
equation, there exists the "superficial" singularity at p = 0. (Of course, there 
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k rho 



Fig.3 The warp factor a'(p)/a, Solution 2 in (j39|) and (^0| 

Horizontal axis: kp. 




k rho 



Fig.4 The warp factor uj'(p)/a, (19). Solution 2 in ( |39|) and 
Horizontal axis: kp. 



Rie alpha^2 




Fig. 5 6D Riemann Scalar Curvature R(p)/a 2 , (||). Solution 2 in (|^) and 

(|4(i|). Horizontal axis: kp. 
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Fig. 6 The Higgs Field P(p)/vq obtained by the Runge-Kutta method. 
Ultra start. Horizontal axis: kp. 
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Fig. 7 The warp factor a'(p)/a obtained by the Runge-Kutta method. 
Ultra start. Horizontal axis: kp. 



is no singularity in the obtained solution.) To avoid the singular point, we 
must start, not from p = but from a point slightly away from p = in 
the positive direction, say, kp = 0.1. In this case, we must properly choose 
the initial values of P(kp = 0.1), a 1 (kp = 0.1) and uj'(kp = 0.1). For these 
values (the vacuum parameters and the initial values) we can borrow the 
information from the analytic result in the previous section. 

□ In Fig.6,7, 

and 8, we plot the results for P(p), cr'(p) and uj'(p) respectively. 

The numerical results well reproduce the solution of the previous analytic 
approach in the ultra-violet region. (Note that we do not assume the form of 
the solution (flS|).) The calculation, however, breaks down at kp ~ 1. This 

14 See App.B (i) Ultra- Violet Start for detail. 
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Fig. 8 The warp factor iv'(p)/a obtained by the Runge-Kutta method. 
Ultra start. Horizontal axis: kp. 



is because the insufficiency of the correct choice of the parameters and the 
initial values. We already know how finely those quantities should be chosen, 
in the 5D case, in order to obtain the required boundary condition |lQ[ . See 
Appendix B for detail, where the results of the infra-start are also presented. 

Note that the fact that the cosmological constant and other parameters 
must be "fine-tuned" to give the required boundary condition, just means 
those are fixed to be some definite values. It is similar to the situation 
the energy eigen values of the hydrogen atom are fixed by the boundary 
condition. It does never mean unstableness of the solution. 
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9 Discussion and conclusion 



We summarize the present result by comparing it with the case of 5D model 
by Randall and Sundrum. 





5 Dim.(Randall-Sundrum) 


o Dim. 


Dimensions of 
Vacuum 
Parameters 


[v ]=(mass) 3/2 
[A]=(mass) _1 
[A] = (mass) 5 


[vo]=(mass) 2 
[A]=(mass)~ 2 
[A] = (mass) 6 


Cond. on 
Cosm. Term 


< a < o 
Lower Bound exists 


A < 
No Lower Bound 


Cond. on 
Fund. Parameters 


7 « i « M 


7«K<M 

Pc 


Asym. Behav. 
of 

Vac. Parameters 
(same k-dep.) 


-A ~ M 3 k' 2 
v ~ M 3 / 2 

A ~ M~ 3 k 2 
a/A^o 2 ~ k 


-A ~ M 4 fc 2 

u ~ M 2 
A ~ M" 4 *; 2 


Mass Relations 


M p i [M 

m ~ y- 

A M = M 2 A 4d ~ -M 3 k 


~ ^ see-saw rel. 
A M = M 2 A 4d ~ — M 4 


4D Cosm. Term 


\J-Aid ~ k ~ M 3 x 10 -38 GeV 


\/-A 4d ~ fc ~ M 2 x KT 19 GeV 



Table 1 Comparison of 5D model and 6D model. 



We add some numerical fact about the see-saw relation (^) . If we take as 
k ~ 10 _41 GeV (inverse of the cosmological size), we obtain M ~ 10~ n Gev= 
10 _2 eV which is the order of the neutrino mass. This is the case 1) in the 
second paragraph of Subsec.5.3, which was abandoned there. Q This choice 
is, however, attractive in that it gives the right value of the cosmological 
constant. If this numerical fact is not accidental and has meaning, it says 
the cosmological size is related to the neutrino mass when it is "see-sawed" 
with the Planck mass. These three fundamental scales are geometrically 

This choice means, as referred in subsec.5.3, the space-time behaves as six dimensional 
at the cosmological scale. In this connection, some interesting discussions are made in 
El E| . They examine some scenario where the world looks like D-dimensions with D>5, 
above a very large (cosmological) size. 
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related. The 6D fundamental scale M gives the mass value, and which 
suggests the neutrino mass is radiatively generated in 6D quantum gravity. 
We consider this is an interesting mechanism to explain the fact that, in 
the log-scale plot of the characteristic mass scales, the neutrino mass is 
located at the middle of the two scales: the Planck scale(10 19 GeV) and the 
cosmological constant scale(10 _41 GeV). Of course, in order to confidently 
accept this view, we should solve the problem presented in Subsec.5.3. 

One of the fundamental parameters, p c , is introduced in Subsec.5.1 as 
the infrared regularization. This is quite natural in the standpoint of the 
discretized approach such as the lattice. The treatment, however, should 
be regarded as an "effective" approach or a "temporary" stage of a right 
treatment. One proposal is made in |l(J for the 5D case where a new infrared 
regularization is introduced. The metric (||) does not have any singularity 
or horizon in the region < p < oo. p = oo can be regarded as a horizon. 
The deficit angle there is 2n. 

We hope the present result of the hierarchy model will lead to the bridge 
between the field theory of the soliton physics and the string field theoy of 
the D-brane physics. 
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10 Appendix A: Calculation of Coefficients 

In this appendix we solve the REE [M = 1) of Sec. 3: 

Eq.(A) 3a" --+a'J + uj" - {Jf + 2- = 2P' 2 ,(12) 
P P 
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Eq.(B) 



8a' 2 + 2 a - 



2a'u' + 2a" = V ,(pt) 



Eq.(C) 4ct" - Wa' z = (P 
in the following form: 
a'(p) = 



J2u 



J{p) = 
P{P) 



oo 

«E 

n=0 
oo 

«E 

n=0 



C2n+1 



(2n + l)! 

^2n.+l 

(2n + 1)! 

e2n 



' 2 -m 2 l- P^ + V) 
P 2 

{tanh(fcp)} 2n+1 
{tanh(A;p)} 2n+1 



m 



v o E ^yy{tanh(M} 2ri ,(19) 



n=0 



where a = V-A/IOM" 2 , @, and V = (A/4)(P 2 



-v 



2^2 



+ A, (P. A, w and 



A are some constants called vacuum parameters in the text. 1/k is another 
parameter called thickness parameter. The integer parameter m in Eq.(C) 
is taken to be m = as explained in (|l^). The coefficient-constants c's, cfs 
and e's are to be determined. 1/p, appearing above, can be expressed as 



1 

P 



2k 



S2r, 



taoh(fcp) ^ (2n)! 



{tanh(A;p)} 



2» 



where 



d 



2n 



In 



1-x . 



S-2,, 



1 

50 = 2 



x=0 

12 



dx 2n 
1 

"3 ' "* 5 

The above expansion formula is a key equation to solve the REE in the form 



S2 



«4 = =" 



(41) 



of (19). All coefficients are finite as shown above. We can take the limit p 
oo above, and see the infinite series converges and gives : 



._ 



(2n)! 



0. 



It is a key, in the present treatment of the infinite series, that we take the 
expansion using powers of tanh p not those of e _2p . 16 

In the following of this appendix, we take the abbreviated notation t = 
tanh(fcp). Using a relation t 1 = tanh(fcp) = k(l — t 2 ), we obtain 



P' = kv (l-t 2 )J2 



e2n 



n=x (2n " 1) 



t 



2n-l 



16 We can not reexpress the RHS of the first equation of ( |4lj ) as 1/p 
with finite coefficients. The same notice is said about (|l9|). 



Eoo 
n = 



(2n)! 1 



-2np 
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oo 



// 7 /-i ,2\ c 2n+l ,2n , r . \^, c 2n+l C 2ra _i u 2 ni 

a =Ml-OL7o^* = fc«{c 1 + ^(-— - )t } , 



(2n)! L1 '^(2n)! (2n-2)! 



re=0 

oo 



7 OO 1 1 

' - Mi - «')E ^ - + - ^f) ■ (42) 

Several useful expansion formulae are listed below. 

/ oo N 

a o; r 1 +2N i i i s 2n C2N-2n+l 

— = 2ka > scUrt where \sc\n = > N , ttttt — ttt 

P ^ J J ^ (2n)!(27V-2n + l)! 

r i 1 r i c 3 s 2 c 3 ci c 5 c 3 ci 

No = -od = 2d. Ni = ^ + ¥ c 1 = ---, [ flC ] 2 = _____, • 

/ oo AT , 

- = 2to^H^t 2iV where = £ " 2n d 2iV-2n+i 



P to ' £-(2n)!(2JV-2n + l)! 

oo AT , 

wV = a 2 [cfc]jv i 2JV+2 where [dc]^ = V / 2ra+ \, , C2Af - 2re+1 - 

^ J L J ^ (2n + l)!(2iV-2n + l)! 

[dc] = dici , [dc]i = -(c?ic 3 + d 3 ci) , • 

OO N 7 7 

c /2 = a 2 V^t 2 ^ 2 where [drf^ = V ,\, ,„ , 

^ J 1 J ^ (2n + l)!(2iV-2n + l)! 

oo AT 

a ' 2 = a 2 £ [cc] ^+2 where [cc] _ y> ggn+l C 2jV - 2ra+1 

jfeo J L J ^ (2n + l)!(2iV-2n + l)! 

[cc] = (ci) 2 , [cc]i = ^ClC 3 , ••• 

oo A" 
j->2 2 r 1 J.2AT i r 1 e2ra e 2N-2n 

P = v £ [ee]* t where [ee]* = £ (2n) , (2jy _ 2n) , 



Ar=0 n=0 

Ho = (e ) 2 , [ee]i=e e 2 , [ee] 2 = ^e e 4 + ^(e 2 ) 2 

oo 

^ 2 = iV(i-t 2 ) 2 Ertt 

A=l 

oo 

= fcVfleeftf 2 + (M' 2 - 2[ee][)t 4 + £ ([«& - 2[ee]'j r _ 1 + M^_ 2 )t 27V } 

AT=3 

AT 

e 2n e 2 Ar_ 2n+2 



where [ee]^ = ^ 



J (2n- 1)! (27V-2n + l)! 
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Mi 



[ee 

N 



2 = ge 2 e 4 



■2\2 



(eoY 



v o }Nt where [e } N = ^ [ee] n [ee]jv-n 

AT=0 ra=0 

1 



[e 4 ]i=2(e ) 3 e 2 , [e 4 ] 2 = -(e ) 3 e 4 + -(e ) 2 (e 2 ) 



6 V ' 2 

Now we find recursion relations between coefficients by inserting expansion 
formulae above into the REE (12,13,10). 

t°-part (c\,di) 



(A) From (12) we obtain 



ka{2ci + 3di} = ,(44) 



4 5 8 5 

ka{-c 3 + -d 3 } - 2fc 2 vo 2 (e 2 ) 2 = fea{-ci + -di} - a 2 (ci - , (45) 

3d 3 3 



t 2 -part (c 3 ,d 3 ,e 2 ) 
t 4 -part (c 5 ,d 5 ,e 4 ) 



7 



13 



1 



-,c\ 



ka < 3 



L 60~ u ' 120 UJ 3' 
+TU d 3 + \d x } - a 2 {l(d lC3 + d 3 ci) - ld!d 3 } - 4k 2 v 2 (e 2 ) 2 , (46) 

t 2N (N > 3)-part (c 2 jv+i, dzN+i, &2n) 

d 2 N-l 



( C 2 N+1 C2N-1 

\(2N)\ (2JV-2)! 



2[sc]iv + 



/ ^2iV+l 

V(22V)! ~ (22V -2)! 



+ 4[sd] 



iV 



+c?{[dc] N -i - [dd]jv_i} - 2A;V{[ee]^ - 2[ee] , JV _ 1 + [ee]^_ 2 } = , (47) 
(B) From ( |i~3| ) we obtain 

i°-part (ci,eo) 



2ka Cl - \\vo\l - (e ) 2 ) 2 = , 
o z 



(48) 



t 2N (N > l)-part (c 2 jv+i,e 2 jv) 



/ca < ( 



.C2JV+1 



C2AT-1 



' (2.2V)! (22V -2)! 
(C) From ([u]) we obtain 



+ 2[sc]7v| - a?{[dc] N -i + 4[cc]7v-i} 
+AV{^ee] 7 v-^[e 4 ]^v} = 0, 



(49) 



i°-part (ci,eo) 
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4fcaci - j\v 4 (l - e 2 ) 2 = A , (50) 
t 2 -part (c 3 ,e 2 ) 

2kac 3 - k 2 v 2 (e 2 ) 2 + -At> 4 e (l - e 2 )e 2 = 4kaa + 10a 2 (ci) 2 , (51) 

i 4 -part (05,64) 
-kacs + ^-^k 2 v 2 e 2 + -^At> 4 e (l - e 2 )| e 4 = 

2kac 3 + ^a 2 cic 3 - 2£; : V(e 2 ) 2 - ^ 4 (1 - 3e 2 )(e 2 ) 2 , (52) 
o o 

t 2Ar (iV > 3)-part (c 2N+1 ,e 2N ) 

ika (S^T ~ (2N N -2)l ) ~ 10a2 t cc ] Ar - 1 - k 2 v 2 {[ee}' N - 2[ee]' JV _ 1 + Hjv_ 2 } 

+iAt !o 4 (2[ee] iV - [e A ] N ) = , (53) 

We have now obtained all necessary recursion relations and are now ready 
for determining all coefficients (c 2n+ i, d 2n+ i, e 2n ). 

i) t°-part, (ci,di,e ) 

Relevant equations are (f5}),( |4q) and (|50|). (^8|) is, however, equivalent to 
(|50|). Hence one coefficient, say eo, remains as a free parameter. 

^SlWd-eo 2 ) 2 + A } , 
oi = --ci , 

eo = — P(iO = 0) : free parameter , (54) 

ii) t 2 -part, (c 3 ,d 3 ,e 2 ) 

Relevant equations are (|45|),(ff9|) with JV = 1 and (|>l]). There are two 
solutions. 



Solution 1 



4, . .a, x 2 x 



e 2 = , c 3 = 2c 1 + 5-(c 1 ) i , d 3 = --(ci + 5-(ci) i ) 



Solution 2 

4 fc 2 



1 A^o 2 0\ 

e2 = --^-e (l - e ) , e 7^ 0, ±1, 
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d 3 = -^(ci + 5|( Cl ) 2 ) . (55) 

Solution 2 includes Solution 1 if we allow eo = 0, ±1. 

iii) t 4 -part, (c 5 ,d 5 ,e 4 ) 
Relevant equations are (p6|),(|49|) with N = 2 and (|52]). They are solved as 

{k 2 v 2 e 2 + ^Xvo 4 e (l - e 2 )}e 4 

4 5 
= ka(--c 3 + 2ci) + a 2 ■ -(dic 3 + d 3 d + 2cic 3 ) 

-\v A ^(l-3eo 2 )e 2 2 + 6k 2 vo 2 e 2 2 , 

5 Aw 4 100 
ac 5 = - — — — e (l - e )e 4 + — — ac 3 + 4aci 
12 k 9 

a 2 10 5 At>o 4 2 \ 2 
j(diC3 + d 3 ci + 8cic 3 ) - -— — (1 - 3e )e 2 , 

120 2 13 

ad 5 = -2ac 5 + —{-kv 2 e 2 e 4 + -^-ac 3 

1 11 , 2 , 

--aci + y^ q «3 + ^aai 



fc 6 

iv) t 2N (N > 3)-part, (c 2N+ i,d 2N+1 ,e 2N ) 
Relevant equations are (|47|) , fl49|) and (p3|), 



1 (dic 3 + ci 3 ci-2d 1 d 3 )-4^o 2 e2 2 } (56) 



, f Q /C2AT+1 C2AT-1 \ . , (d 2N+ i d 2 N-l \ . ,r ,i 

to \ 3 1 w " (2F^2)lJ - 2[sc]n + { W¥ " W^yJ + 4[s(% 

+a 2 {H7v-i - [dd] N -i} - 2k 2 v 2 {[ee]' N - 2[ee] , w _ 1 + [ee]' N _ 2 } = . 
ka \^^NS\ ~ (27^- 2)! ^ + 2 t sc ] iV } ~ « 2 {[ dc ]^-i + ±[<An-i. 

+\v o *{l[ee} N -l[e%} = 0., 

4ka ~ (27^- 2)! ) ~ 1 ° Q;2 [ cc ] Ar - 1 _ fc 2 «o 2 {[ee]jv - 2[ee]' Ar _ 1 + [ee]^_ 2 ] 

+^o 4 (2[ee] 7 v - [e A ] N ) = . 
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Because [sc]jv, [sd]jv, [ee]jy, [ee]^ and [e 4 ]N have the terms C2JV+1, cfejv+i and 
C2N linearly (for N > 3) and other parts of the relevant equations are made 
of lower-order terms, they are three coupled linear equations for the three 
variables C2N+1 > <h,N+l and eiN- Hence we conclude that (C2JV+1, cfoiv+i; e 2iv) 
are fixed by the lower-order coefficients. 



11 Appendix B: Numerical Analysis by Runge- 
Kutta 



We explain the numerical analysis of Sec. 7 of the text. Instead of the func- 
tions a',u' and P, it is convenient to use normalized ones: 



a a 



They satisfy the unit boundary condition: P, £, Q — > 1 as p 



(58) 

00. In terms 



of these functions, we can rewrite the REE equations (12, 13,10) as 
-S = --=-E + ^(Sfi + 4E' Z ) + — (p2 



ld„ 1„ . ,„, v lAu 4 M- 4 ,~ Q 1AM" 4 



fc tip 



, + f (££7 + 4£ 2 ) + 
kp k 



8 fca 



d 2 + 



k ap kp kp k 8 



2 fca 

x9 1AM" 4 

1) H ; 

; 2 fca 



+ 



aM 4 1 a 2 M 4 



k 2 VQ 2 



4 Jfe 2 



1) 2 + 



A 



k 2 VQ 2 
\ 1/2 



This is the standard form of the Runge-Kutta equation for P, £ and SI 
with the coordinate kp. All dimensionless constants appearing in the RHSs 
are rewritten, in terms of the dimensionless vacuum parameters in 
yf^K/M 2 k = A, \v 4 /M 4 k 2 = L, v 2 /M 4 = B, as 



-AM" 



ka 



'10A 
aM 4 1 A 



Xv 4 M- 



ka 



kvtf 



10^ 



W 
k 2 



10- 

L 
B 



L a _ A 

A ' k = 7To 

-A A]_ 
' fc 2 V ~ B 
a 2 M 4 _ 1_^2 
/c 2 u 2 ~ 10 5 



(60) 



29 



The values of (A, B, L) must be properly obtained corresponding to a given 
boundary condition. For the case P(kp = 0) = eo = —0.8 case (Solution 2), 
we have obtained their values as (A, B, L) ~ (ao, A)> <5o) = 
(1.864786,0.6424929,15.90691) in Sec. 6. We take these values in this ap- 
pendix. As for the initial condition for (P, X,f2), we analyze two cases. 

(i) Ultra- Violet Start 

First we numerically solve the differential equation in the direction from 
the small kp value to the large one. Because the origin kp = is the 
"superficial" singularity in (|59|), we must take, as the initial point, some 
point slightly away (in the positive direction) from the origin. Say, kp = 0.1. 
As the boundary value of P(0.1), X(0.1) and $7(0.1), we take the values 
obtained from the analytic result of Sec. 6. They are -0.791057,-0.123938 
and 0.082700, respectively. The results are shown in Fig.6,7,and 8 of Sec. 7. 
Those figures surely reproduce the analytical result of Sec. 6 up to kp ~ 1. 
It strongly supports the correctness of the solution found in Sec. 3 and 6. 
It also means the validity of the truncation approximation to determine 
the vacuum parameters. The calculation stops with producing a imaginary 
value. It occurs at the place where the inside of the square root in the last 
equation of (|5^) approaches but is still required to be non-negative. The 
present precision is insufficient to keep non- negative in such nearly region. 

(ii) Infra-Red Start 

Next we numerically solve the equations in the direction from the large kp 
value to the small one. Say, from kp = 3.1 to kp = 0.1. As for the initial 
values P(3.1), E(3.1) and 0(3.1), one idea is to take the analytic result of 
Sec. 6, as in the (i) case. It, however, fails. The more stringent requirement 
of the positivity, than the case of the UV start case, is demanded. We can 
see the situation by writing the inside of the square root of ([59]) as 

-—p) 2 = - JL^±e + + 5n)(E - n) 

kdp J y/VQBkp 5 B K A ' 

-^(l-^) + I|(P 2 -l) 2 >0 . (61) 

In the asymptotic region kp — > do, we know P — > 1, S — ► 1, Q — * 1, H/kp — > 
0. Therefore the above quantity becomes quite small and, at the same time, 
is required to be non-negative. The failure means the approximation of the 
3-terms truncation used in Sec. 6 is insufficient in the IR region. Here we 
cannot help being content with the qualitative correctness of the solution. 
In Fig. 9, 10, 11, we show the results for the initial data: the initial point 
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Fig. 9 The Higgs Field P(p)/vq obtained by the Runge-Kutta method. 
Infra start. Horizontal axis: kp. 
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Fig. 10 The warp factor a\p)ja obtained by the Runge-Kutta method. 
Infra start. Horizontal axis: kp. 



kp = 3.1; the initial values P(3.1) = 1.0, £(3.1) = 1.12,0(3.1) = 1.12. They 
are consistent with the analytical result in the qualitative way. 
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